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ABSTRACT 

Gamma-ray bursts are believed to be produced when fast moving ejecta from 
some central source collides with slower moving, but relativistic, shells that were 
ejected at an earlier time. In this so called internal shock scenario we expect 
some fraction of the energy of the burst to be carried by slow moving shells 
that were ejected at late times. These slow shells collide with faster moving 
outer shells when the outer shells have slowed down as a result of sweeping up 
material from the ISM. This gives rise to a forward shock that moves into the 
outer shell producing a bump in the afterglow light curve of amplitude roughly 
proportional to the ratio of the energy in the inner and the outer shells. In 
addition, a reverse shock propagates in the inner shell and produces emission at 
a characteristic frequency that is typically much smaller than the peak of the 
emission from the outer shell by a factor of ~ 77q^(-E2/-Ei)^'^5 and the observed 
flux at this frequency from the reverse shock is larger compared to the flux from 
the outer shell by a factor of ~ 8(7oc-E'2/-E'i)^^^; where 7oc is the bulk Lorentz 
factor of the outer shell at the time of collision, and Ei & E2 are the total 
energy in the outer and the inner shells respectively. The Lorentz factor is 
related to the observer time as ~ 5{t/dayY^^. The shell collision could produce 
initial temporal variability in the early afterglow signal. The lack of significant 
deviation from a power-law decline of the optical afterglow from half a dozen 
bursts suggests that E2/E1 is small. Future mult i- wavelength observations 
should be able to either detect bumps in the light curve corresponding to both 
the forward and the reverse shocks or further constrain the late time release of 
energy in ejecta with small Lorentz factor, which is expected generically in the 
internal shock models for the gamma-ray bursts. 



Subject headings: gamma rays: bursts - relativistic shock 
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1. Introduction 

The recent mult i- wavelength observations of GRB afterglow (Costa et al. 1997, van 
Paradijs et al. 1997, Bond 1997, Frail et al. 1997) support the fireball shock model for 
GRBs (Wijers et al. 1997, Vietri 1997, Katz & Piran 1997, Waxman 1997). According to 
this model GRBs are produced as a result of internal shocks when a fast moving shell runs 
into a slower moving shell that was ejected at an earlier time. The relative kinetic energy of 
motion of the shells is converted into observed gamma-ray emission via a relativistic shock. 
The observed afterglow emission in this scenario is produced when the shell slows down 
as a result of interaction with the ISM. This internal-external fireball model (Meszaros & 
Rees 1997, Sari & Piran 1997a-b) requires a complicated central engine, which operates 
for a long time (as long as the duration of the burst) producing a highly variable flow. 
The observed temporal structure seen in GRBs, follows, to a large extend, the temporal 
structure produced by the source (Kabayashi, Piran & Sari, 1997). Meszaros et al. (1998) 
considered the effect of an inhomogeneous fareball, where the fireball parameters had 
angular dependence, on the afterglow. At the end of the internal shock phase we are left 
with a rather ordered flow in which faster (merged) shells are the outermost ones and slower 
shells follow behind them (see fig. 1). If the source continues to operate for even longer 
time, producing even slower shells, they will follow behind the earlier ejected faster moving 
shells. 

As the fast outer shells move outwards they begin to interact with the ISM and 
decelerate. Eventually the slower inner shells will catch up with the outer shells. Thus one 
should expect late interaction of slow shells with faster shells that has been slowed down. 
Our goal in this paper is to examine such collisions and to provide some observational 
signatures of such interactions which should be seen in future afterglow observations. The 
current and the future observations can also be used to place a limit on the fraction of 
energy coming out in slow moving ejecta at late times. As an example the hypernova model 
(Paczynski, 1997), in which perhaps most of the energy comes out from the central engine 
in the from of slow moving ejecta that lags behind the outer fast material, should show 
deviation at late times from smooth power-law decline of the after-glow light curve. 

Rees & Meszaros (1998), Panaitescu et al. (1998), and Cohen & Piran (1999) have 
recently considered the continuous limit of the delayed emission problem focusing on the 
forward shock emission. They discussed a source that emits a wind whose characteristics 
varies as a power law with time. However, in view of the internal shock model we expect 
that an impulsive situation in which two shells collide should be more relevant to GRB 
afterglows. We consider therefore, a toy model consisting of two shells that undergo 
collision, and we calculate synchrontron emission from both the forward as well as the 
reverse shock that propagates into outer and the inner shells respectively. The interaction 
of the outer shell with the ISM is described well by the adiabatic Blandford-McKee (1976; 
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hereafter BM) self-similar solution. The inner shell catches up with the outer shell when the 
outer shell slows down and its Lorentz factor becomes approximately equal to the Lorentz 
factor of the inner shell. The shell collision produces two new shock waves: a forward shock 
that moves into the outer shell and a reverse shock that propagates into the inner shell. We 
analyze these shocks and estimate the resulting synchrotron emission. 

If the collision takes place not too early during the afterglow the system is adiabatic 
(this is probably a valid approximation in GRB afterglow from about half an hour after the 
burst; Granot, Piran & Sari, 1998) and the outer material has arranged itself with a BM 
self-similar profile. In this case the overall effect of the collision between the outer shell, 
with energy Ei and the inner shell, with energy E2 would be a transition from one BM 
solution, with a total energy Ei to another one, with a total energy r^ Ei + E2, i.e. the 
observed flux at a fixed frequency will increase by a factor of (1 -|- E2/ EiY''^. This estimate 
is not exact as there would be some enhanced energy losses form the shocks that arises 
during the collision. At most wavelengths we should expect a smooth transition from one 
solution to the other. The exception is at those frequencies at which there is significant 
emission from the reverse shock. As we show later the emission from the forward shock is 
rather similar to the emission from the shock that arises from the interaction of the outer 
shell with the ISM. The total emission from the reverse shock (the shock that propagates 
into the inner shell) is smaller but typically comes out at a much lower frequency. At these 
frequencies the impact of this additional emission is likely to be very pronounced. 

Early on, during the early afterglow we expect that the system is radiative (it must be 
radiative during the GRB phase, otherwise the energy budget would be much higher). In 
this case we expect that the outer shells that has collided with the ISM would cool rapidly 
and the collision between the inner and the outer shell would be between two cold materials. 



The calculations presented in section |2.1| and ^]2| are valid only within the adiabatic regime, 
namely several hours after the beginning of the afterglow. However, similar qualitative 
behavior also applies during the early afterglow while the shock is still radiative. In 
particular, after the internal shocks have arranged the outflow in a monotonically increasing 
function of radius, the successive collisions of these shells would lead to a variability in the 
early afterglow with the variability time-scale of the order of the time since the explosion. 

Note that we have made several simplifying assumptions in our model. In addition 
to the above mentioned assumptions of local spherical behaviour (which is a good 
approximation at this early stage even if the system is beamed) and adiabaticity (which 
might be broken as we discuss below) we have assumed that (i) the ratio of electron and 
magnetic field energy density to the total thermal energy is the same behind all shocks and 
(ii) that the ratio of specific enthalpy in the inner and the outer shells is approximately 
given by the ratio of their energy. This latter assumption is valid only if the ratio of the 
thickness of the slow shell and its radius is approximately given by the inverse of the square 
of the Lorentz factor. We discuss these assumptions after equation ^ in section ^ and in a 
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footnote in section pJ] . Clearly our results should be modified if these assumptions are not 
satisfied. 

We describe our physical model in section 2. A detailed discussion of the shock 
conditions, when a cold and a hot shells collide, is presented in §2.1, and in §2.2 we discuss 
the effect of the stratification of the outer shell, using the Blandford-McKee (BM) solution, 
on the shock propagation and emission. We discuss briefiy the situation in a radiative case 
in §2.3 The synchrotron emission from the shocks and its observational consequences are 
discussed in §3. Implications to recent observations and predictions for future observations 
are contained in the final section. 



2. The Physical Model 

Our basic model consists of a cold inner shell that is colliding with an outer hot shell. 
The outer shell is slowing down as a result of collision with the interstellar medium (ISM) 
so that the slower moving inner shell eventually catches up with it. 

Let us take the energy in the outer shell, as measured in rest frame of the center of the 
explosion, prior to the shell collision, to be Ei, and the energy in the inner shell to be E2. 
The Lorentz factor of the outer shell, 70 (t), decreases with time. At late times when 70 (t) 
is much less than the initial value, 70J, the Lorentz factor is given by 



where po is the density of the interstellar medium, and 

1/3 



3^1 



(2) 



is the time when 70 (t) = 7oi/2. 



The radius of the inner shell at time t is Riit) ~ t — 1/(2^1), and the radius of the 
decelerating outer shell is -Ro(^) ~ t — t/[8'jQ{t)]; 74 is the Lorentz factor of the inner shell, 
and all quantities are measured in the rest frame of the center of the explosion. The inner 
shell runs into the outer shell when R^i^t) = Ro{t), and at this time 70 (t) = 70c = 74/2, 
and thus the relative Lorentz factor of collision is approximately 1.25 so long as we neglect 
the slowdown of the inner shell. The time when the shells collide is tc ~ ^0 (701/74) ^'^^- ^ot 
the self-similar structure of the outer shell, described by the Blandford-McKee solution, 
it can be shown that the Lorentz factor of the inner shell w.r.t. the rest frame of fiuid of 
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the outer shell, when the collision takes place, is about 1.25. Therefore, the relative speed 
of propagation of the two shocks, one propagating in the outer and the other in the inner 
shell, is time independent and is equal to 0.6 c. In this case the outgoing forward shock, 
that propagates into the already shocked material of the outer shell, is weak and the reverse 
shock propagating into the cold inner shell is strong but mildly relativistic. 

The space can be divided into five different regions. (I) The ISM is the outermost 
region which consists of cold gas, is taken to be at rest relative to distant observers, and is 
characterized by a single parameter the density, Uq. (II) The outer shell and the shocked 
ISM consists of relativistic electrons and protons moving with a bulk Lorentz factor of 
7o(t) relative to the ISM; the thermodynamical variables in this region are denoted by a 
subscript 1 and the equation of state is pi = ei/3. Region III consists of material of the 
outer shell that has been heated by the forward shock resulting from shell collision; this 
region is also relativistic and the variables are denoted by a subscript 2. (IV) Part of the 
inner shell that has been heated by the reverse shock; electrons are relativistic in this region 
but the protons are not and consequently the equation of state is ps = 4e3/9. (V) This 
region contains unshocked cold inner shell characterized by density n^. This shell moves 
with a Lorentz factor 74 w.r.t. the ISM and 714 relative to region 2. The different regimes 
are depicted in Fig. 2. 



2.1. Shock conditions 

We analyze a pair of forward shock and reverse shock formed when a cold region 
(^4777,^0^ ^ 64) collides with a hot relativistic region {nirripC^ ^ ei = 3pi). The system is 
characterized by three parameters: the energy density in the hot region, ei, the particle 
density in the cold region, n^, and the relative Lorentz factor between the inner and the 
outer shells, 714. In fact most of the quantities are determined just by the dimensionless 
ratio of the enthalpy density w^jwi = 3n4mpC^/(4ei) ~ E2/E1, and by 714. [| A fourth 
parameter (ni) does not appear in the shock conditions but it determines the density in the 



^In deriving the relation w^/wi ~ E2/E1 we made use of the fact that the thickness of the 
cold inner shell is Ra/jI as seen in an inertial frame. The thickness of a shell, in shell rest 
frame, which is undergoing radial expansion at the sound speed ~ c/S^^^, is approximately 
ctg = ct/'y^; where tg and t are the time since the shell was expelled from the central source in 
the shell reast frame and as seen by an inertial observer respectively. Thus the shell thickness 
in an inertial frame is ~ ct/74 provided that the shell radius at the time of collision is much 
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shocked region III (^12). 

Two shocks form, the forward shock propagates into the hot outer shell and the reverse 
shock propagates into the cold inner shell. Each of the shocks satisfy three jump conditions 
which ensure the conservation of particle number, energy and momentum across the shock. 
For the forward shock propagating into region II (see fig. 2) we have: 

ni'jivi = n2'j2V2 (3) 

Wi'jfvi = W2llv'2 (4) 

Wxl\v\ + Pi = ^272^2 + Vi (5) 

where 71 (72) and the corresponding velocity fi (^2) are measured with respect to the rest 
frame of the shock and n^, e,, p, and Wi are the particle density, thermal energy density, 
pressure and enthalpy measured with respect to the rest frame of the fluid in region 'i+1'. 

The forward shock propagates into hot material that has been already shocked when 
the first shell interacted with the ISM. The mass energy density in the outer shell is 
negligible relative to the thermal energy density i.e. riimpC^ <ti ei, and the equation of state 
is: pi = ei/3. As the shocked material can be only hotter, similar conditions should hold 
in region III i.e. n2mpC^ <C 62 and p2 = 62/8. We also expect that 62 > Ci. With these 
conditions we find: 

"2 = :^ (6) 



-,, = 2%±£2) p) 



6ei 



72 



/3(ei + 3e2 



862 \ (87? - 9) 



9(7? 



_ (3ei + 62) (ei + 3e2) _ 3.? - 1 
^^ V 16eie2 2v^ ^ ' 



n2 = rii-fivl 



It - 9) 



^\ 7i-l \ei(3ei + e2 



'' ^'' + ^^^).ni (10) 



larger than its radius at expulsion, and at collision the sound speed of the inner shell is not 
highly sub-relativistic. The former condition is satisfied at all times during and after the 
7-ray burst. The shell thickness is less than ct/74 when the latter condition is violated and 
in this case the peak frequency of emission from the reverse shock is correspondingly smaller. 



7t2 



62 

n2mpc'^ 



In 



(87? - 9)7? 
\i 9(7? -1) 



;ii) 



where 7/ is the relative Lorentz factor between regions II and III, thus it is the shocked 
matter velocity in the rest frame of the unshocked material, and 7^2 is the "thermal" 
Lorentz factor of protons in region III or the typical Lorentz factor of the random motion 
of the protons and is crucial in determining the emission from this region. 

The reverse shock is different. In the region of the reverse shock (marked IV in 
fig. 2) the electrons are relativistic but the protons are non-relativistic or at best mildly 
relativistic. The equation of state in this region is ps = fje^/S, where 77 = 1 for relativistic 
protons and 77 = 4/3 for non-relativistic protons. For the unshocked cold shell 64 = p4 = 0. 
Separating regions III and IV is a surface of contact discontinuity. Both regions move 
with the same velocity and both have the same pressure: p2 = Ps- The energy densities 
in the two regions are related by 62 = rje^. Using these relations and the jump conditions 
analogous to equations (3)-(5) we obtain: 



74 
73 
Ir 

63 

7t3 



(627/ + 7l3 - 714) (e2 + 8714) 
\ (62 (-1 + 3r]) + 8713 - 8714) 7Z4 



(62 + na) (62 + 723 - ^4) 
^ (I862 + 9723) (562 + 97Z3 - 9714 



(7762 + n^) (62 + 8714) 

^ [(1 + 877) 62 + 8713] 7l4 



(1 + 67/) 77,4 + y 47/627Z4 (1 + Si]) + (1 + 6r]) ni^ 

2 

^62 

(63/7^3772^6^). 



(12) 
(18) 
(14) 

(15) 

(16) 
(17) 



Here, 7^ is the Lorentz factor of the motion of the shocked matter in region IV relative to 
the cold matter in region V. To avoid cumbersome equations we have expressed 74, 73 and 
7j. in terms of 773. However, 773 is given in a closed form in terms of 62, and 774. So 62 is the 
only remaining unknown parameter in equations (7)-(17). Since regions HI and IV move 
with the same velocity we can express the relative Lorentz factor between the cold inner 
shell, region V, and the outer shocked material, II, in terms of the Lorentz factors 7/ and 



714 = 7r7/ + V'[(7.'-1)(7|-1)] (18) 

For given parameters (61, 774 and 714) we can solve equations (9), (14), (15) & (18) to 
determine 62 and thus all other variables in the shocked regions. Note that 62/61 depends 
only on 714 and the ratio 774/61 oc 7^4/^1. 
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The ratio of thermal energy densities in regions III and II, 62/61, is a function of the 
relative Lorentz factor of collision (714) and w^/wi. For a fixed value of 714 ~ 1.25 we find 
from the numerical solution of the conservation equations that 62/6I ?» [1 + [w4^/wi)^'^^]. 
The ratio of the densities n2/ni is an even weaker function of [104^/ Wi). 

There is a region in the parameter space for which the above set of equations have no 
solution. The no-solution case physically corresponds to when the enthalpy density, or the 
energy, in the inner shell is too small to drive a forward shock into the outer shell. This is 
best seen in fig. 3 which depicts 62/61 vs. W4/W1 for several values of 714. Since regions 
III is shocked, 62 > 61 with 62 = 61 a limiting case in which the forward shock disappears. 
Therefore, for a given value of 714 the minimal value of W4/W1 is obtained when 62 = 61. 
The minimal value of W4/W1 for a forward shock to occur increases with decreasing 714; for 
7i4 = 1.25 the minimum value of w^/wi is 0.36. 

The correct treatment of the problem in those cases where we don't have a solution of 
the conservation equations is to relax the homogeneous shell approximation used thus far 
and use the Blandford-McKee solution for the structure of the outer shell. As ei, the energy 
density in region II decreases with distance from the interface of regions I & II, there will 
always be a point in region II at which w^/wi is large enough so that forward shock forms. 
This is discussed in the next section. 



2.2. Structure of the outer shell and shock solution 

We turn now to determine the parameters for the outer shell (region II of fig. 2). This 
region contains shocked material from the ISM and the baryonic ejecta from the explosion. 
There is a shock between this region and the cold ISM. In a simple relativistic shock 
between two cold shells we have (see e.g. Piran 1999): 

ni ^ 47ono, 61 = AjQUompC^. (19) 

The simplest approximation would be to assume that region II is homogeneous and to 
use equation (19), together with the adiabaticity condition: 

E, = {An/3)R''nompC^^l (20) 

to determine the conditions in region II. However, the inner shell overtakes the outer one 
when the outer shell has shocked on the ISM and has slowed down to a Lorentz factor that 
is roughly equal to that of the inner shell. At this stage a rarefaction wave has propagated 
through the outer shell and it has settled down to a Blandford-McKee (BM) self-similar 
solution where the density, enthalpy, bulk Lorentz factor etc. decrease with distance from 



-10- 



the outer surface of the shell. The collision of the shells and the resulting structure and 
emission from the shocked regions should be calculated using the stratified structure of the 
outer shell given by the BM solution. 

The BM solution is expressed in terms of the similarity variable x, which is defined by: 

X = 1 + 167? (^) , (21) 

where R is the radius of the shock front, r is the distance inward from the shock from 
(measured in the ISM rest frame) and 7j is the Lorentz factor of the matter just behind the 
shock: 

RVf = ^/^ (22) 

where /^ = E/[{47!'/3)nQmpc'^]. R and 7/ are related to the observed time at which radiation 
from the front of the shock reaches the observer as: 

^"'^ ^ 68 /^ "^ Ta^ ^^^^ 

The BM self-similar profile is given by: 

n' = 47^noX-'/' , 7 = T/X^'^' , e' = 4nom,c'^jx~''/'\ (24) 

where n' and e' are the number density and the energy density in the local rest frame, 
respectively, and 7' is the Lorentz factor of the bulk motion of the matter at radius (R — r). 

A complete and exact solution of the colliding shell problem requires determining 
the propagating of shock wave through a BM stratified medium. This could be done 
numerically. However, an approximate analytic solution should be sufficient for our purpose. 

As discussed at the beginning of §2 the relative Lorentz factor of collision i.e., the 
relative Lorentz factor of region V and the bottom of region II (714), is 1.25 and this is 
independent of x- Thus, as the forward shock propagates outward through the stratified 
outer shell, the value of 714 remains nearly constant. However, the value of W4/W1 decreases 
as x^'^^^'^j W4/W1 at X = 1 is approximately time independent. 

In the last sub-section we found that for a fixed value of 714 the forward shock dies out 
if W4/W1 falls below some critical value; for 714 = 1.25 the minimum value of w^/wi is 0.36. 
Thus, if the energy of the inner shell {E2) is less than about 0.36 Ei, the forward shock 
stalls and turns into a rarefaction wave when it reaches the surface Xc ^ 0.5{Ei/E2y^^^'^ ■ 

For Xc ^ 1, the forward shock traverses through the outer shell and the energy density 
at a radius Xj within the shell, increases as a result of the shock by a factor proportional to 
X^^^^^, and the particle number density increase is an even weaker function of x- Since the 



-11- 



synchrotron emission scales as e^/n, most of the emission from the outer shell, pre- as well 
as post-collision, is generated near the top of the shell. The overall increase in the emission 
from the outer shell can, therefore, be calculated by treating the shell to be homogeneous 
with values of ni, ei and 71 corresponding to x = 1 in the BM solution. Emission from the 
reverse shock can also be calculated using the same values of these parameters. 

For Xc ^ 1 the forward shock is very weak, it stalls at x = Xci and has small effect 
on the emission from the outer shell. The emission from the reverse shock can however be 
still much larger than the emission from the outer shell at the frequency corresponding to 
the peak of the synchrotron radiation from region IV. This emission can be calculated by 
applying the results of the last section with 714 = 1.25, W4/W1 = 0.36, and other parameters 
in region II corresponding to x = Xc of the BM solution. 



2.3. Radiative Collisions 

Slower shells that catch up with faster decelerating shells during the first few hours of 
the burst find the outer shell to be cold as the synchrotron cooling time is shorter than the 
dynamical time, and thus we need to consider collision between two cold shells. 

The calculation of shocks resulting from the collision of two cold shells resembles the 
calculation of energy emitted from internal shocks (Sari & Piran, 1997a). It is also rather 
similar to the calculation given in §2.1. The main difference is that the pressure in region 
II is zero, and as a result there is always a forward and a backward shock. At early times, 
within a few hours of the burst, the shock is radiative and the Lorentz factor drops off with 
time as t~^ . The relative Lorentz factor of the shells at the time of collision is about 1.5. 
Simple kinematic calculation suggest that, equal masses colliding with a relative Lorentz 
factor of 1.5 and radiating efficiently will radiate away ~10% of their total energy. The 
rest of the energy will simply be added to the kinetic energy of the ejecta plowing through 
the ISM. We find that the Lorentz factor of the forward and the reverse shocks are about 
1.12 when ni ~ n4. Both shocks are similar and they emit their energy at a much lower 
frequency than the outermost shock that is propagating into the ISM at this stage. We 
expect, therefore, that this emission will be a very significant contribution to the emission 
at this lower wavelength at this stage. 



3. Synchrotron Emission 

The peak of synchrotron emission occurs at frequency Vm ~ QellB / {mec) , and the 
frequency integrated emissivity is given by e si o'rcne7g-B^/87r; where Qe-, fUe, 7e are electron 
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charge, mass, and thermal Lorentz factor respectively, B is the magnetic field in the fiuid 
rest frame, Ug is the electron number density, and ax is the Thomson cross-section. The 
emission at low frequencies (z/ <^ //„) drops off as u^^^ and at high frequencies {u ^ //„) 
scales as z/~(*'~^)/^, where p Ri 2.4 is the power-law index for the energy distribution of 
electrons. 

If the energy in the magnetic field is taken to be some fraction, ^b, of the thermal 
energy density then fi^/Svr = ^B^e'^ef^e- And thus the ratio of the peak synchrotron 
frequencies in regions III and II is 



Y2rn _ fe2\ 



5/2 



ni 



2 nU-5 



H - 9) 



(25) 
_n2. " 

and the ratio of the observed flux from the outer shell after the forward shock has traversed 
through it, and the flux in the absence of collision, at a frequency much greater than the 
peak of the emission is given by 

hi:^) =!,[-) [—) . (26) 

where 7j = 7j(l + f/fo) is a factor by which the Lorentz factor of the outer shell increases 
as a result of shell collision. The dependence of the observed flux on 7j can be more rapid 
than the linear function considered above, depending on the temporal proflle of the shell 
acceleration. 

At late times, when the collision has run its course, the shells have merged and settled 
back to a BM solution, the observed flux is proportional to {Ei + E2Y'^^^''^'^] in the absence 
of the shell collision the flux would have been proportional to Ef . Thus the increase 
in the observed emission due to the forward shock is approximately (1 + E2/ EiY'^^'^'^f^. 
Numerical calculation of flux increase, using equation (26), yields result that is consistent 
with this estimate. 

The reverse shock traversing into the inner shell is mildly relativistic but very strong. 
The density enhancement in this case is about 4, i.e. n^/ui ~ 4. Regions IV and III are 
separated by a surface of contact discontinuity across which the pressure is continuous 
but the density is not. The thermal Lorentz factor of electrons in region IV, 743, can be 
calculated from the continuity of pressure at this interface and is given by 



_ 63 _ 3mp 
rrien^ 40me 



Wie2 



^461 



(27) 



and the thermal Lorentz factor of electrons in region II, 7^1, can be shown to be equal 
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to (7oc/2^/^)(mp/me) provided that electrons and protons are in equipartitionj^ where 
7oc = 7o(^c) is the Lorentz factor of the outer shell at the time of collision. The ratio of the 
peak synchrotron frequencies in regions IV and II is given by: 



Vim 27/1 Vei/ 77oc Vei/ Kw^ 

We have made use of the relation B3/B1 ^ (3e2/4ei)^''^ in deriving the above equation; Bi 
& i?3 are magnetic fields in regions II & IV respectively. We have furthermore assumed that 
the fractional energy density in the magnetic field and electrons in regions IV and II are not 
different. This is an economical assumption which eliminates having to introduce additional 
free parameters. Moreover, there is some observational support for this assumption from 
the analysis of GRB970508 (Weijers & Galma, 1998) and the prompt optical flash of 
GRB990123 (sari & Piran, 1999). It is straightforward to modify the results presented here 
if this assumption is found to be invalid. 

Using the result of §2.1 we can write 62/61 = f{wi/wi) ~ [1 + {wi/wiY''^^]. Moreover, 
it can be easily shown that w^/wi ~ E2/E1. Making use of all these results we find that 
the Lorentz factor of electrons in region IV is smaller than the thermal Lorenz factor 
of electrons in the outer shell by a factor of /3 ~ 57oc {E2/E1) f^^{E2/Ei). And so the 
peak frequency of the synchrotron emission emanating from the shocked inner shell is 
smaller than the characteristic emission frequency from the outer shell by a factor of 
P'f-'^' - Hc{E2/Eif\ 

The ratio of observed flux from the inner shell, after the reverse shock has propagated 
through it, and the outer shell in the absence of collision, at the peak frequency 
corresponding to the inner shell, is given by: 

hjl^Zp) _ 7/^3^3 \Vlrn] ^''^ _ ifH f 537*1 ] ^^^ /29) 



7/^3 

ni 


537tl 
.5i7i3_ 



/i(«^3p) niBi Lz/3„ 

Since rii = 2^^'^'jocnism, Wi ~ 87Q^njsm'^p/3, "^,3 ^ An^ ~ 4w4/mp, we find that 
ns/rii ~ jocWi/wi. Substituting this into the above equation we find that the observed 
emission from the inner shell at frequency 700^^3, is larger than the emission from the outer 
shell by a factor of ~ 8['yocE2/E- 



15/3 



As an example consider that the inner shell is ejected with a Lorentz factor of 5, and 
its energy is comparable to the energy in the outer shell then the radiation from the shocked 



^The equipartition assumption is in fact not needed since we are only interested in the 
ratio of peak frequencies and emissions from regions II, III, and IV. 
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inner shell comes out at a frequency which is smaller compared to the peak frequency 
of emission from the outer shell by a factor of almost 10^ and the observed flux at this 
frequency is dominated by the inner shell by a factor of about 10^; if the energy in the inner 
shell were 1/5 of the energy in the outer shell then the frequency ratio would be about 30 
and the flux from the shocked inner shell larger by a factor of about 8. 

Fig. 4 shows the peak frequency and the emission from the reverse shock, and fig. 5 
shows the modification to the light curve due to the forward shock. 

The synchrotron cooling time in shell rest frame is ts = Sircme/icTTB^le), and the 
dynamical time of the shell is td = t/70. Thus the ratio of the cooling to the dynamical 
time for the outer shell is given by: 



ts 87rcme7o ml 1.58x10 H^f^^ 



2 1 c;cin-4+i/2 

(30) 



td (TrB^-fet ieiBCTTcmlriQilt ^^^buJ^eH^' 

where ^e & ^b ai'e the fraction of the total energy density in electrons and the magnetic field. 
For ,^e ~ 'Cb = 0-1) ^0 = 1 and E^ji = 1 the outer shell becomes adiabatic approximately one 
hours after the burst in the observer's frame. The synchrotron cooling time for the inner 
shell heated by the reverse shock is even longer because the magnetic field strength in the 
regions III and IV are approximately equal as a result of pressure balance across the surface 
of contact discontinuity (see fig. 1), and 73^ ^ 724 ~ 70- Thus our estimate of energy flux 
calculated above under the assumption of adiabatic shock is valid as long as we restrict 
ourselves to shell collision after about one hour. 



4. Discussion and comparison with observations 

In the scenario that gamma-ray bursts are produced by internal shocks it is expected 
that some fraction of the energy of the explosion is carried by ejecta that is moving with 
moderate Lorentz factor which does not collide with faster moving material that was ejected 
at earlier time until the faster shells have been slowed down by the ISM. We have explored 
the consequences of this possibility in this paper. The goal is to provide some constraint on 
the temporal behavior of the energy release from the explosion. 

When a slower moving shell hits a decelerating faster shell from behind, it results in a 
correlated increase of emission at all frequencies. The relative Lorentz factor of collision can 
be shown to be about 1.25, and the amplitude of the enhanced emission depends on the 
ratio of the energy in the inner and the outer shells as well as the observed frequency. 

The enhancement to the observed emission from the forward shock, that propagates 
into the outer shell, is approximately (1 + E2/ EiY''^ and the characteristic synchrotron 
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frequency is slightly larger than the frequency of the outer shell in the absence of shell 
collision. 

The emission from the reverse shock, which propagates into the inner shell, is at a 
frequency that is smaller than the peak frequency of emission from the outer shell by a 
factor of 77Q(tc)(-E2/-Ei)^'^; where 7o(tc) is the Lorentz factor of the outer shell at the time 
of collision. The emission from the inner shell, at the peak frequency, is larger than the 



emission from the outer shell at the same frequency by a factor of ~ 8{'yQE2/E 



|5/3 



Thus the observed spectrum, after shell collision, at a fixed time is expected to have 
two peaks. The lower frequency peak arises from the inner shell and the higher frequency 
peak from the outer shell. Moreover, the light curve at a fixed frequency has a bump that 
can be very dramatic at low frequencies. The results presented in this paper are subject to 
certain assumptions described in the introduction viz. the fraction of the thermal energy 
carried by electrons and magnetic field is the same in all shocks, and the thickness of the 
inner shell is of order its radius divided by the square of its Lorentz factor which is expected 
from causality considerations. 

The observed light curve of the gamma-ray burst GRB970228 appears to have a 
plateau during 6-10 days after the burst with an amplitude of about 0.5 mag in the R-band 
(Fruchter et al. 1998). If this were to arise as a result of shell collision, one would conclude 
that the energy in the inner shell is about 40% of the outer shell. 

The light curve of the x-ray afterglow of GRB970508 shows a complicated behavior. 
At about 4 days after the burst the light curve shows some deviation from a power-law 
fall-off (Sokolov et al. 1997), from which we infer the late time energy injection to be 
less than about 10% of the initial burst. At early times the behavior was more dramatic. 
The observed x-ray flux in the 2-10 kev band fell between 11 and 16 hrs after the initial 
gamma-ray burst and then rose sharply by almost a factor of two from about 16 to 20 hours. 
After a gap in the observation from 20 hours to 60 hours this was followed by a powerlaw 
decline with time as ~ t~^ (Piro et al. 1999; |astrp-ph/9902013D . A possible interpretation 



of the sharp rise is that the outer decelerating shell, which is producing the x-ray afterglow 
resulting from the shocked ISM, is being hit from behind by a shell with energy roughly 
equal to the energy of the outer shell. If this is the correct interpretation then we expect 
that the optical emission will follow this overall trend. After an initial decline of flux that 
was seen in the early optical observation (the earliest of which was about 5 hours after the 
burst) a sharp increase in the optical flux began one day after the GRB cf. Sokolov et al. 
(1997). The optical flux peaked around two day after the burst and from then on it began a 
powerlaw decline. It seems that there is a lag of almost a day between the optical and the 
X-ray light curves, which is not expected in the shell collision model as both the X-ray and 
the optical emissions are expected to arise from the forward shock (the emission from the 
reverse shock should peak in the milli-meter wavelength at this time, and at this wavelength 
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the emission from the reverse shock is larger than the emission from the forward shock by 
about two orders of magnitude). However the gaps in the X-ray and optical observations 
do allow a much smaller lag or even no lag at all. If so than the jumps in the observed 
X-ray and optical fluxes are consistant with enhanced emission due to a colliding shells 
(Panaitescu, Meszaros & Rees, 1998) with a significant energy release. 

The detection of optical flash associated with GRB-990123 with peak magnitude of 
8.95 in the V-band 50 second after the initial gamma-ray burst (Akerlof et al. 1999) was 
remarkable. The recent work of Sari and Piran (1999) provides a nice fit to the optical and 
radio observations based on this model. It is interesting to note that these observations can 
also be explained as a result of collision of two shells which collide with a moderate relative 
Lorentz factor of order 2 or less. Such collisions are expected in internal shock models when 
either a somewhat faster moving shell takes over a slower shell or when a slower moving 
shell catches up with a faster moving but decelerating shell. In the early, radiative, phase 
of the shock the Lorentz factor of the outer shell decreases with time as 1/t^, and in this 
case it can be easily shown that when slower moving shells catch up with a faster shell 
their relative Lorentz factor is 1.5. We have analyzed a simple model consisting of wind 
from the central source with randomly fluctuating Lorentz factor, lasting for about 100 sec 
(the duration of the GRB), and find that the total fluence in the optical wavelength band 
during the gamma-ray burst is about 0.1% of the x-ray and the gamma-ray fluence, which 
is consistent with the observation reported by Akerlof et al. (1999). 

The smoothness of the observed light curves of other GRB after-glow emission suggests 
that the late time release of energy, in slower shells, is typically quite small. 

Future, multi-wave length, observations should provide more stringent condition on the 
ejection of moderate/low Lorentz factor material from the central engine of the gamma-ray 
bursts at late times. 
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Figure Captions 

Figure 1. — The panel to the left shows the Lorentz factor of ejected shells as a function of 
time (random function). The panel to the right shows the Lorentz-factor-distribution as a 
function of radial distance from the center at some time when the slower moving shells have 
been taken over by faster shells leading to a monotonically increasing Lorentz factor with 
distance (a constant was subtracted from the distance). The width of hues is proportional 
to shell mass. 

Figure 2. — Shows the various shock regions when a cold shell collides with a hot shell. 

Figure 3. — The graph show the ratio of the thermal energy density in regions (III) and II 
(^2/61) as a function of the ratio of the enthalpy density of the cold inner shell, region V, 
and the outer shell, region II, (^4/^1); the ratio w^/wi is approximately equal to the ratio 
of energy in the inner and the outer shells. The four curves corresponds to four different 
values of the relative Lorentz factor of the two shells (714) : the value of 714 for the bottom 
to the top curves are 1.25, 1.5, 2.0 and 3.0 respectively. Note that the minimum value of 
W4/W1 for which there is a forward shock i.e. 62/61 > 1, decreases with increasing 714. The 
value of 7i4 is close to 1.25 for the case where a slow moving shell catches with a faster but 
decelerating shell. 

Figure 4. — The top panel shows the ratio of the peak synchrotron frequency in the region 
of the reverse shock to the outer shell multiplied by 7oc (70c is the thermal Lorentz factor 
of gas in the outer shell or region II). The lower panel shows the ratio of emission from 
the reverse shock region and the outer shell, at the characteristic synchrotron frequency of 
region IV, divided by 7q^ , as a function of W4/W1 ~ E2/E1 at a time when the reverse 
shock has reached the back end of the inner shell. The solid curve is for 714 = 1.25 and the 
dotted curve is for 714 = 2.0. The horizontal part of the curve for 714 = 1.25 corresponds 
to the forward shock stalling at a point in the outer shell where w^/wi = 0.36, and the 
physical parameters in the reverse shock is determined by this critical point. 

Figure 5. — The effect of shell collision on the afterglow light curve (schematic). The 
continuous curve corresponds to E2/E1 = 0.4 and the dashed curve is for E2/E1 = 0.7 



